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ABSTRACT 
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Major Professor: Steven A. Stage, Ph*D* 

The effects of external parameters on the surface 
heat and vapor fluxes into the marine atmospheric boundary 
layer (MABL) -during cold-air outbreaks are investigated 
using the numerical model of Stage and Businger (1981a). 
These fluxes are nondimensional ized using the horizontal 
heat (gi) and vapor (g a ) transfer coefficient method first 
suggested by Chou and Atlas (1982) and further formulated 
by Stage (1983a). 

In order to simplify the problem, the boundary layer 
is assumed to be well mixed and horizontally homogeneous, 
and to have linear shoreline soundings of equivalent 
potential temperature and mixing ratio. Modifications of 
initial surface flux estimates, time step limitation, and 
termination conditions are made to the MABL model to 
obtain accurate computations. 
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The dependence of g x and g a in the cloud topped 
boundary layer on the external parameters (wind speed, 
divergence, sea surface temperature, radiative sky 
temperature, cloud top radiation cooling, and initial 
shoreline soundings of temperature, and mixing ratio) is 
studied by a sensitivity analysis, which shows that the 
uncertainties of horizontal transfer coefficients caused 
by changes in the parameters are reasonably small. 
Therefore, the surface heat and vapor fluxes are also 
weakly dependent on the external parameters, and can be 
wel 1 -estimated by the horizontal transfer coefficient 
method even from rather poor measurements. 


It 

is 

found numerically in 

this paper that 

the 

concept 

of 

horizontal transfer 

coefficients can 

be 


applied to both cloud-free and cloud-topped regions of the 
marine boundary layer for the estimation of surface fluxes 
at any given fetch. Although it is best to estimate 
fluxes by direct computation using the numerical MABL 
model whenever possible, values of g x and g a are close to 
those obtained analytically for the dry case. Thus, the 
analytical solutions of the horizontal transfer 
coefficients can be used to obtain surface flux estimates 
which are accurate enough for many purposes without 
actually numerically integrating the MABL model. 
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CHAPTER 1 INTRODUCTION AND MOTIVATION 


Cold-air outbreaks often occur after a frontal 
passage during winter since large temperature gradients 
exist in the frontal area* Once the cold front passes over 
warm water, the large temperature and vapor contrasts 
between the warm water surface and the cold and dry arctic 
air result in abnormally strong surface heat and vapor 
fluxes into the boundary layer* These fluxes are critical 
factors in the boundary layer modification and are 
important as energy sources for atmospheric motions above 
the boundary layer* Atlas et al* (1983) suggested that 
the large energy input to the atmosphere during cold air 
outbreaks associated with mesoscale convergence zones, and 
shoreline shape might be important in cyclogenesis* 
Numerous studies show that brief cold air outbreaks 
supply a large portion of air-sea fluxes* Hence, they are 
very important in determining the energy budgets of sea 
and atmosphere during winter, late fall, or early spring* 
Chou and Atlas (1982) first suggested that for any 
given fetch the net sensible and latent heat fluxes per 
unit travel are proportional to the surface temperature 
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and mixing ratio differences between sea and shoreline. 

Based on the assumptions of well mixed boundary 
layer and horizontal homogeneity in the shoreline 
direction. Stage (1983a) derived analytic solutions for 
the fluxes before the cloud edge. His solutions provide 
a nondimensional method to parameterize the surface heat 
and vapor fluxes from warm water surfaces of dry 
(cloud-free) boundary layer during cold air outbreaks. 

The goal of this paper is to extend Stage's 
nondimensional method (1983a) of computing dry boundary 
layer surface fluxes into the cloud-topped marine boundary 
layers during cold air outbreaks. A modified numerical 
version of the boundary layer model, which was first 
developed by Stage and Businger ( 1981a, b) for cold air 
outbreaks, is presented in this paper and used to compute 
the nondimensional flux parameters, horizontal heat (gx) 
and vapor (g 3 ) transfer coefficients for both cloud-free 
and cloud-topped regions. Sensitivity analysis of how gi 
and g a are affected by the uncertainties of external 
parameters measured by satellite and shoreline sounding 
is the major approach of this study. 

From the results of this work, a better understanding 
of the changes of heat and vapor fluxes during the layer 
evolution is obtained. Most of all, it is found that the 
uncertainties of g x and g a caused by the changes of 
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external parameters are not very large. Therefore, the 
surface heat and vapor fluxes at any given fetch can be 
reasonably estimated by the horizontal transfer 
coefficient method, for most of the cold air outbreak 
episodes even with very poor measurement. 

The technology of today can measure some of the 
boundary layer parameters at any point of the fetch by 
satellite, as discussed by Allison (1984). These 
parameters include surface wind speed or surface stress, 
integrated cloud top temperature, water surface 
temperature, and integrated liquid water and water vapor 
content in the atmosphere. By using the satellite 
measured data and a well-mixed boundary layer assumption, 
Allison's method (1984) can calculate the surface heat, 
vapor and momentum fluxes at any point of the fetch. 
However, the work done in this thesis provides a 
complement to Allison's work. Given shoreline sounding 
measurements and similar boundary layer assumptions, the 
surface heat and vapor fluxes between the shore and any 
given point can be computed. Thus, the technique verified 
in this thesis provides an alternate way to compute 
fluxes, which can be used when shoreline soundings are 
available, but not the satellite parameters needed for 


Allison's method 
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Chapter 2 introduces the basic thermodynamics and 
equations of the Stage and Businger's model ( 1981a, b); and 
the horizontal transfer coefficient method proposed by 
Stage (1983a) for the parameterization of heat and vapor 
fluxes to the cloud edge during cold air outbreaks is also 
presented. gi and g a computed by the modified numerical 
marine boundary layer model are compared with the analytic 
solutions of horizontal transfer coefficient method. The 
result shows that the numerical and analytical solutions 
are identical to each other in the dry (cloud-free) 
boundary layer. 

In chapter 3, the modified numerical model is applied 
to estimate the boundary layer evolution into the cloud- 
topped regions. The sensitivity analysis shows that 
changes of heat and vapor fluxes relative to the 
fundamental boundary layer parameters are not very large. 
This suggests that the horizontal transfer coefficient 
method is a promising way of estimating surface fluxes at 
any given fetch, regardless of the large measurement 
uncertainty of satellite and sounding data. 

Finally, in chapter 4, the conclusions of this paper 
and outlook for further research work will be presented. 



CHAPTER 2 FORMULATION OF THE MODEL 


2.1 Introduction of The Model Physical Background 

In this chapter, the model developed by Stage and 
Businger ( 1981a, b) to simulate the well-mixed, cloud- 
topped planetary boundary layer over warm water during 
cold-air outbreaks will be discussed. Such outbreaks 
freguently occur in many parts of the world during the 
fall and winter, especially over the western sections of 
midlatitude oceans off the east coasts of Asia and North 
America. In these regions, cold-air outbreaks combined 
with large air-sea temperature differences, high winds, 
synoptic scale convergence, and relatively cold dry air 
can cause much greater heat and moisture fluxes than 
oceanic means. So far, other studies of similar boundary 
layer situations include analysis of AMTEX data of cold 
air masses which leave China and cross the warm water of 
the Kuroshio, (Ninomiya 1974, 1975; Ninomiya and Akiyama 
1976); Chou and Atlas's (1982) studies of cold-air 
outbreaks over the warm water area, off the south shore 
of Long Island during the winter of 1979; and Stage and 
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Businger's (1981a, b) studies of cold-air outbreaks over 
Lake Ontario from IFY6L data. 

The model mentioned above for simulating the 
evolution of the cloud-topped marine boundary layer during 
cold air outbreaks is further modified and applied in my 
work. This model predicts all of the mean thermodynamics 
of the boundary layer, such as the rate of change of mean 
eguivalent potential temperature, mean mixing ratio, and 
boundary layer depth following a column of air as it 
traverses a warm water surface . 

All these changes in mixed-layer properties result 
from the following effects: surface turbulent fluxes of 
heat and water vapor, radiative warming of cloud bases and 
cooling of cloud tops, condensation and evaporation, 
entrainment of dry warm air from the inversion base, and 
lifting (subsidence) imposed by convergence (divergence) 
in the mixed layer. These processes are shown 
schematically in Fig. 2.1.1. 

Parameterization of net sensible and latent heat 
fluxes per unit travel from the surface, and their 
evolution during cold air outbreaks are the main concern 
of studies in this paper. This parameterization theory 
was first developed by Chou and Atlas (1982). They 
noticed that the net surface sensible heat flux per unit 
travel , H„/x, at any given distance from the shore , 
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fetch, is approximately proportional to the surface 
temperature difference between the sea and shore. Also 
the net surface latent heat flux per unit travel , H L /x, 
at any given fetch is approximately proportional to the 
surface mixing ratio difference between the sea and shore. 
They then suggested that these surface flux per unit 
travel terms can be simply expressed as horizontal 
transfer coefficients multiplied by the surface 
temperature or the mixing ratio differences. 

Later in this paper, more details about how the 
horizontal heat and vapor transfer coefficients are 
affected by the fundamental parameters of the atmospheric 
boundary layer will be described. 


2.2 Fundamental Variables and Basic Thermodynamics 

The model is based on the assumption that the 
equivalent potential temperature $ m and total water mixing 
ratio q T are linear functions of height at the shore and 
are well mixed throughout the boundary layer. Therefore, 
the layer can be treated as a slab. Sharp jumps exist in 
the profiles of $ m and q T at the top of the mixed layer 
because of the inversion. Also, it is assumed that the 
boundary layer is horizontally homogeneous in the 
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shoreline direction. 

The vertical profiles of $ m and q T , and the 
corresponding vertical fluxes W'q T ' are shown in 

Fig. 2. 2.1. All the variables used in this model are 
defined below. 

The independent variable is either: 

t time, 

or 

x fetch. 

The fundamental dependent variables are: 

Z. ensemble mean mixed layer depth due to 

entrainment; 

Z 0 ensemble mean mixed layer depth due to 

convergence; 

0 m vertical and ensemble mean mixed layer 

equivalent potential temperature; 

q T vertical and ensemble mean mixed layer- 

total water vapor mixing ratio. 

At the shoreline , = o, and Z. = Z D + Z. = Z. , where Z B 

is the ensembled mixed layer depth, and the initial values 
are specified for Z., 8 and qT* Then all other 


10 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Fig. 2. 2.1 Profiles of idealized MABL model. (From 


Stage and Businger, 1981a, Fig. 3) 

a. Equivalent potential temperature. (^.) 

b. Total water vapor mixing ratio, (q-r) 

c. Equivalent potential temperature flux. 

d. Total water vapor flux. <W'q-r‘) 
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variables of the boundary layer can be calculated from 
these fundamental variables. If applied to data, ensemble 
averages must be estimated by time averages in the steady- 
state condition; averages parallel to the shoreline with 
homogeneous upwind condition in that direction; or areal 
averages with small extent in the streamwise direction. 

The basic thermodynamic properties of the model 
fundamental variables are shown below. The equivalent 
potential temperature is defined as : 

6 m = 0exp( Lq„ / C*T ), (2.2.1) 

where L is the latent heat of evaporation of water; 

Cp is the specific heat of air at constant 
pressure; 

q v is the water vapor mixing ratio. 

An accurate approximation of (2.2.1) can be written as: 

6rn = 9 + ( L/C* ) q v . (2.2.2) 

The total water mixing ratio is defined as: 


q T = 


qv + qx , 


(2.2.3) 
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where 

q v is water vapor mixing ratio, 
qi is liquid water mixing ratio. 

The Clausius-Clapey ron equation may be expressed as: 


de = ( 6eL/RT„ a ) dT rt , (2.2.4) 

where 

6 is the ratio of the molecular weight of water to 
dry air ( € S 0.622 ), 

Td is the dew point of the air, 

e is the water vapor partial pressure of the air 
which can be expressed as: 

e * Pq v /6 , (2.2.5) 

where P is pressure of the air. 

By differentiating (2.2.5) we have: 


dq„ / q„ 


de/e - dP/P, 


( 2 . 2 . 6 ) 
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and by the use of (2.2.4), the ideal gas law, and the 
hydrostatic equation , the following equation is derived : 

dq„/q v = ( 6 L ) / ( RTa a ) dT* + (g/RT)dZ. (2.2.7) 

He can integrate (2.2.7). If the saturation mixing ratio 
(q r ) corresponding to the dew point temperature (T dr ) is 
known, then the pair of related functions are inferred 
as : 

Q(T ) = q r exp ( -i^(JL-~) ), (2.2.8) 

and 

Tdo(q) = ( + "t^T >"** (2.2.9) 

where 

Q(T) is the saturation water vapor mixing ratio 
for air at Z = 0 and temperature T; 

Tdo(q) is the dewpoint for air at Z = 0, and water 
vapor mixing ratio q. 

Notice that these two functions are inverses of each 
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other, thus 

q = Q (Tdo ( q ) ) , 


(2.2.10a) 


and 


T = T„ 0 (Q(T)). (2.2.10b) 

For air below its saturation level, q„ is conserved, so 
that dq v = 0, and from (2.2.7) we can derive: 

dTd = - (gT/L6)dZ (2.2.11) 

by using T = T*. 

Integrating (2.2.11) as a parcel is lifted from the 
surface up, and assuming T = constant, 

T„ = Tdio ( q v ) - (gT/€L)Z. (2.2.12) 

From (2.2*12), we see that for nonsaturated conditions 
the approximate dew point lapse rate is -(gT/6L) = 

1.76 C/km. We can define dew point potential temperature 
as : 


= 


Te* + 7* Z 


(2.2.13) 
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By combining (2.2.12) and (2.2.13), we have: 

9* = Tae(q v ) + ( t - gT/feL)Z. (2.2.14) 

Another useful quantity is liquid water potential 
temperature, 0i, which is defined as: 

0i = 6m - (L/CsOq,. (2.2.15) 

For adiabatic parcel motions 0i is conserved because 9m 
and q T are conserved. So, from ( 2 . 2 . 2 ) , ( 2 • 2 . 3 ) and 
(2.2.15), we know: 

0i = 9. for qi = 0 (2.2.16) 

(i.e. for nonsaturated air parcel, qi=0 ). 

The level at which dewpoint and temperature become 
equal is defined as the lifting condensation level, Z e * 
This can be derived from ( 2 . 2 . 13 ) , ( 2 . 2 . 14 ) , ( 2 . 2 . 16 ) and 
0i = T + /'Zc by assuming the parcel has q x = 0 at level 
Zc* Thus, 


L dO 


(gT/6L)Z c 


0i ” /^Zc , 


(2.2.17) 
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Ze = ( 0i-T* o )/( Y - gT/6L ). (2.2.18) 

The virtual potential temperature, , is defined as: 

= 0 ( 1+ (l/6-l)q v + qi) (2.2.19) 

= ( 0 . - (L/C*)q v )( 1 + (1/6 - 1 )q„ + qi). 

This temperature includes the effects of both water vapor 
buoyancy and liquid water drag. 

The last important thermodynamic property to be 
introduced is the dewpoint virtual potential temperature 
defined as: 

= *«.( 1 + (1 /6 - l)q v + qi). (2.2.20) 


2.3 Analytic Profiles of Heat and Vapor Fluxes in the Dry 
MABL Model 

If divergence is negligibly small then the net 
sensible and latent heat input from a warm water surface 
to the atmosphere in cloud-free regions (commonly referred 
to as the dry boundary layer), can be computed by vertical 
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Integration from the surface to the inversion layer of 
the change of and q T in soundings. This computation 
method is adequately only when the atmospheric divergence 
is negligible, because there is hardly any change of the 
sounding profiles above Z B . Thus, we can express the net 
sensible and latent heat input to the boundary layer as : 



where 


(2.3.1) 


(2.3.2) 


H is the net sensible heat input. 

Hi, is the net latent heat input. 

The analytic solution for a dry boundary layer was 
first developed by Stage (1983a). The more complicated 
condition of a cloud-topped boundary layer (called as wet 
boundary layer) after cloud formation during boundary 
layer evolution has not been solved analytically. In this 
section, the analytic profiles of heat and vapor fluxes 
before the cloud formation point will be presented. 
First, it is known that there is no liquid water (qi=0) in 
this region, hence q T = q„. Therefore, the virtual 
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potential temperature and virtual heat flux input to the 
boundary layer can be written as: 


6 ~ = 0( 1 + ( 1/6 - 1 )g T ) 


(2.3.3) 


Let 


H v = SCM - 0„| )dZ . 

K *«° 


(2.3.4 ) 


Substituting (2.3.1), (2.3.2) into (2.3.3) gives 


H v - H + ( 1/6 - 1 )9o(Cp/L)H l . 


(2.3.5) 


Next, neglecting molecular diffusion and horizontal 
convergence of turbulent heat fluxes, the heat 
conservation equation may be written as: 


0 . 


(2.3.6) 


Integrating (2.3.6) with respect to Z, and applying 
boundary conditions leads to: 


“dT 


-~r (W0„'lo + ), 

■^8 


(2.3.7) 


where the boundary condition for the dry boundary layer at 
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Z = Za- and Z = 0 gives: 


( = W * fl v * | O 


at Z = 0, \ 


W*^’ < 


>(2.3.8) 


V. = H'dv'l ZB . = -W.a0 at Z = Zb-/ 


where 


W. is the entrainment rate, 

aQ is the temperature jump at the cloud top. 

Now, take the z -derivative of (2.3.6). Since we assume 

the layer is well mixed, ffv = Q m , and = 0 is a good 

approximation for this boundary layer model. Thus, 
(2.3.6) becomes: 

= Const. (2.3.9) 

Ball (1960) proposed a relationship between cloud top 
entrainment rate and surface heat flux: 




A W '0 V ' I o . 


(2.3.10) 
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A 2 0.2 is a widely accepted value of the entrainment 

coefficient . 

From (2.3.8) and (2.3.9), we know the heat flux of 
boundary layer: 

= VT?7*'Io(1-Z/Zb) + W'0 v 'L»-(Z/Z b ) (2.3.11) 

Further discussion of nondimensionalization of the 

net sensible heat and vapor fluxes per unit travel, and 

the analytic solutions of the evolution of these fluxes 

from real sounding data will be presented in the following 

sections of this chapter. 

For the wet boundary layer case (i.e. after cloud 
formation), the heat and vapor fluxes become much more 
complicated. The radiation cooling rates (R a ) at the 
cloud top, and the radiation warming rates (R c ) at the 

cloud base have to be considered. In order to simplify 

the problem , R a and R c are assumed to be concentrated in 

a very thin layer near the cloud top and base and entirely 

within the mixed layer. Thus, the effect of the thickness 
of the radiative flux divergence layer near cloud top and 
base can be neglected (Stage and Businger, 1981b). 
Through a proof similar to the one above, by using the 
diffusion equation for the total water mixing ratio (qx) 
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and the equivalent potential temperature ( $ m ) r we can 
derive : 


dq T _ . 

dt *z “ °> 


(2.3.12) 


d6« __ aw^e’ e(j£) __ 

~dl dz * dz ~~ 


(2.3.13) 


where 


»(&> 

dZ 


= 0 


everywhere except at Z c and Z 8 , 


and 


iim 

5->o 


JEl 

ft? 


Zc+ J 
Z e -J 


= -Rc 


(2.3.14a) 


Om. 


Fa 2b+j 



(2.3.14b) 


By taking the z-derivative of (2.3.12) and (2.3.13) to 
show that the turbulent fluxes must be linear in the cloud 


and 

subcloud 

regions, and 

then twice 

integrating 

the 

result using 

the boundary 

conditions , 

(2.3.14a) 

and 

(2.3 

. 14b ) , the 

fluxes of total water mixing ratio (W'q T ') 

and 

equivalent potential 

temperature 

(W’0.‘) can 

be 


expressed as (Stage and Businger, 1981b): 
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W * q T * = ( 1-Z/Z. ) Wq T '| 0 - ( Z/Z. )W.Aq T , (2.3.15) 


( = (1-Z/Z. )We’Jo-Z/Z.(-W.A0+R.-R c ), \ 


W 'Q m ' 


for 0 < Z < Z 


C, 




(2.3.16) 


(1-Z/Z.) (W’0. , |o+Rc)+Z/Z.(-W.a^+R»), 


for Z c < Z < Z., 


These Idealized profiles of W'q T ' and W '0.' are shown in 
Fig. 2.2.1 (c), (d) of the previous section. 

In the cloud layer all the air motions are assumed to 
be saturated. Taking the perturbations of (2.2.2), 
(2.2.3), (2.2.19), q„ ' = (dQ.AT/dT )$ ' , and neglecting the 

small terms . 

(2.3.17) 


This gives: 


W'0 V • = B W • Qm ' - 8 W ‘ q T * 


for Zc< Z < Z 


(2.3.18) 
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where 


t 


s ( 1 + — -^5 w 1 + JL.-43&.) -i 

€ dt Hl q>^F' 


has a typical value of 0.5 to 0.6. 


Combining (2.2.19), (2.3.15), (2.3.16) and (2.3.18) gives: 


/=(l--S)W , ^ v ’|o + «J( R .-Rc)-^A i W., 


W’0 V 


for 0 ^ Z < Z c , 


(2.3.19) 


= (1- -j) (f (W^ + R c )-0Wq^) + ^Ra--|AaW., 


for Z c ^ Z < Z m. 


where 

Ai = &B i + ( 1/6 - 1 ) <3t/ 

Aa = ^ 8 - - ( 1/6 - 1 ) fl^qt. 


The idealized profiles of 6>* and W'fl v ' are shown in 
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Fig. 2.3.1a and b. 

Further discussions of the parameterization of the 
net sensible and vapor fluxes per unit travel for the dry 
boundary layer (i.e. with no clouds) will be presented in 
the next section. 

In this paper, the major concern is to use a 
numerical method which can simulate the more complicated 
and realistic evolution of the net sensible and latent 
heat fluxes from a warm water surface into cloud-topped 
region during the boundary layer development. Then, the 
numerical results will be compared to the analytical 
solutions to show the credibility of the numerical method 
used in this paper. 


2.4 Flux Transfer Parameterization in the MABL 

A parameterization method for the heat and vapor 
fluxes of the dry boundary layer, no cloud case, will be 
discussed in this section. This method was first proposed 
and solved analytically by Stage (1983a). For simplicity, 
this method assumes linear shoreline soundings. They are: 


e 


= $m 1 + "TffmZ 


( 2 . 4 . 1 ) 
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and 


*3t ,b* = *3x1 + T^t^b • 


(2.4.2) 


In the last section, we have mentioned the rate of 
changes of surface sensible heat and vapor fluxes. The 
time integrated sensible and latent heat fluxes during 
boundary layer evolution can be written as functions of 
the fetch as follows: 


H = CPCWU) 


Hi, = (PL/U) 


\ ^qT'l, 

J o 


dx 


dx 


H„ = (fCp/Q) 


J a 


dx 


where 


(2.4.3) 


(J is the mean wind speed in the boundary layer, 

x is the fetch of boundary layer evolution as 
measured in the direction of U. 


Chou and Atlas (1982) noticed that for any given 
fetch the mean net sensible heating per unit travel (H v /x) 
is proportional to the surface temperature difference 
between air over land and sea (9„o-9 v i), and that the mean 
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net latent heating (Hl./x) per unit travel is proportional 
to the surface mixing ratio difference between air over 
land and sea (g 0 *qi ) • The analysis plotted by Chou and 
Atlas (1982) is reproduced in Fig. 2. 4.1. 

Stage (1983a) parameterized H „/x and H b /x by using 
the horizontal transfer coefficient terms proposed by Chou 
and Atlas (1982). The net sensible and latent heating 
per unit travel are expressed as the horizontal transfer 
coefficients multiplied by (0vo-5vi) and (go-9i)* 

Hw/x = C t ( 0 VO - ) gi / (2.4.4) 

Ht/x * C a ( qo - qi ) 9a / (2.4.5) 


where 


Ci, C a are the constants for sensible and latent 
heat fluxes respectively. 

The parameters of eg. (2.4.4) and (2.4.5), Ci, C a , gi and 
g a , have been derived by Stage (1983a) as shown below: 

Ca = C^C Ti , (2.4.6) 


C a 


= -^PLC T i , 


(2.4.7) 
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where 

Cn is the value of C T corresponding to an air 
temperature of 9*, = and a sea surface 

temperature of 0 vO . i.e. the surface temperature 
difference between air at shore and sea. 

Cs> is the specific heat of air at constant pressure, 

L is the latent heat release of water vapor 
evaporation. 

The horizontal transfer coefficients of sensible heat, g x , 
and latent heat, g a , are 


9i = X , (2.4.8) 

g a = ( 2 ( 1 + A ) (l-F/ty) ( l -\) +\)\ a /X . (2.4.9) 


The nondimensional terms included in the above equations 
are defined as follows: 

(A) >1 : nondimensional mixed layer height, - which is 


defined as: 
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7l = Z/Z»«, 


(2.4.10) 


where Z» is a scaling height, i.e. the depth of 
boundary layer when fetch approaches infinity. 
It is defined as: 


Z m = ( ‘ ) ( O “ 1 ) /Tv r 

I +4 


(2.4.11) 


where 

o : virtual potential temperature of the sea 

surface, 

9vi : virtual potential temperature of the 

shoreline surface, 

■p v : virtual potential temperature lapse 

rate in the boundary layer, 

A : entrainment coefficient ( A - 0.2 ). 


(B) X s nondimensional slope of mixing ratio, g T , 
profile, which is defined as: 

X = (r./r > v)((&'0-&.i)/(qo-q»)) - (2.4.12) 


where 
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JT, is the nixing ratio lapse rate in the BL, 
go is the mixing ratio at the sea surface, 
qi is the mixing ratio at the shoreline surface. 

(C) X : nondimensional fetch of boundary layer 

evolution, which is defined as: 


X = x/x» , 


(2.4.13) 


where x» is a scaling length in the fetch 
direction. It is defined as: 


Xw = Zboo/ ( ( 1 + A ) Ct x ) , 


(2.4.14a) 


from (2.4.11) we have 


I + 3A , 8vo — 8 vi N 
x - = TvC T i 5 


(2.4.14b) 


(D) F : nondimensional function. 


F 


n.-( 


- ( 1 - y + 


ty-y^) 

A 



(2.4.15) 


? 


where y = 1 - Y£ 
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For 7l« 1 (i.e. Z << Zsoo), during the initial stage of 
the MABL evolution, (2.4.15) contains small differences 
between terms and is hard to compute. However, we can 
apply the binomial expansion in the limit situation 
( 71 << 1) to simplify the problem, such that: 

F s (• • (2.4.16) 

By substituting (2.4.16) into (2.4.9), we can prove that 
9, = #* = gi, for 71 << 1. This has been shown in the 
numerical results for the initial stage run of this paper. 
The related plots will be presented in the later sections 
of this chapter. 

From the definitions of horizontal sensible heat 
transfer coefficient g* and horizontal latent heat 
transfer coefficient g a , (2.4.8), (2.4.9), and (2.4.13), 
we can see that gi and g a are functions of fetch (x). 

How g t and g a change with fetch under different 
conditions is the major topic of the research work in this 
paper. Further details will be discussed in the next 
chapter. 

From (2.4.8) to (2.4.15), nondimensional ization of 
vertical sensible and latent heat fluxes in MABL, all the 
parameters that closely affect the horizontal transfer 
coefficients for the dry boundary layer case are: 
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0 vO ! sea surface virtual potential temperature. It 
is computed from sea surface temperature f0 o ) 
and mixing ratio (g T o)* 

0vi : surface virtual potential temperature at shore. 

It is computed from shoreline surface 

temperature ( Sx ) and mixing ratio (g T i). 

P„ : lapse rate of virtual potential temperature. 

It is computed from the lapse rate of 

temperature ( P 5 ) and mixing ratio ( Hit ) • 

Di : the divergence within the boundary layer. 

D a : the synoptic scale divergence. 

and 

U : the wind speed. 

All of the above parameters have been assumed constant to 
simplify the situation. From (2.3.6) to (2.3.19), and 
(2.4.3) to (2.4.5), we can see that once clouds form, gi, 
g 2 are also affected by: 

R b : the radiation cooling on cloud top. 

R a * — ( Tzb - 4 “ T b k v 4 ) , 

-r C P 

where T B kv is the radiative sky temperature. 

R c : the radiation warming at cloud base. 

Rc * J?Cp ( * - Tzc -. 4 ), 

where To is the sea surface temperature. 

In conclusion, the method of deriving g x and g a with 
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nondimensional parameters provides a means of evaluating 
heat and vapor fluxes with much less information than is 
required for the complete analysis* This method will be 
called the "Horizontal Transfer Coefficient Method”, and 
the solutions of g x and g a can be solved either 
analytically or numerically. 

As learned in this section, the advantage of the 
analytical solutions is that simple observations at shore 
and sea are adequate to solve for g* and g a at any given 
fetch. However, the biggest disadvantage is that this 
solution provides the values of heat and vapor fluxes 
before cloud formation only. 

A schematic illustration of the analytical and 
numerical methods of solving the surface heat and vapor 
fluxes is shown in Fig. 2. 4. 3, where g x and g a are 
analytically solved if the surface temperature, mixing 
ratio over land and ocean, and their lapse rates at shore 
are known. The analytical solutions of g x and g a are 
universal for all the cloud-free boundary layer cases. 
Thus, the surface heat and vapor fluxes per unit travel 
can be derived from (2.4.4) and (2.4.5). It is also seen 
from Fig. 2. 4. 3, g* and g a are numerically solved by 
applying the MABL model (Stage and Businger, 1981). There 
are 9 external parameters needed to run the simulating 
boundary layer evolution model, it is good for both the 
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cloud-free and cloud-topped boundary layer. For the dry 
boundary layer, the numerical solutions of gi and g a are 
identical to the analytical solutions, which will be 
proven in the next section. For the wet boundary layer, 
the numerical values of g A and g a are considered to be 
reasonably accurate. Hence, the surface fluxes per unit 
travel from shore all the way to cloud-topped regions can 
always be estimated by running the numerical MABL model. 

2.5 Comparisons Between The Numerical and The Analytical 

Solutions of The Model 

Stage (1983a) analytically solved g x and g a of MABL 
evolution during cold air outbreaks for the dry (cloud- 
free) case. In order to compute the heat and vapor fluxes 
for both the wet (with cloud) and dry (no cloud) cases, a 
moidified numerical method from Stage and Businger's MABL 
model is applied in this section. This method will also 
be evaluated and compared with the analytical method later 
in this section. 

The numerical method uses the same assumptions of 
linear conditions as the analytical model, given earlier 
by (2.4.1) and (2.4.2). Specification of the initial 
values of the fundamental variables enables the 
calculation of all the variables of the boundary layer. 
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Then the boundary layer evolution is determined by 
integrating the numerical model at every time step. These 
numerical solutions agree well with the analytic solutions 
in the cloud free region. They also provide relatively 
accurate forecasts of the heat and vapor fluxes for the 
cloud-topped region as will be shown in the next chapter. 
A comparison between the analytic and numerical solutions 
of horizontal transfer coefficients is presented in this 
section to examine how well the numerical model works. 

The sounding being used here for the comparison of 
the two methods is the New York sounding (Stage 1979/ 
1981b/ 1983b; Chou and Atlas 1982), taken on February 
1979. It is considered to be a typical air modification 
situation during cold air outbreak episodes. The details 
of this sounding are shown in Chapter 3, Table 3.1.1. A 
vertical profile of the fundamental parameters of the 
boundary layer for the New York sounding is shown in Fig. 
2.5.1 (Stage, 1979). 

Although the results of the two methods shown in 
this section are from one particular sounding, it will be 
proven in Chapter 3 that the horizontal transfer 
coefficients determined by this numerical model are not 
very sensitive to changes in sounding parameters. 
Therefore, some general conclusions can be made with 


regard to the model results based on 


single sounding 



Fig. 2. 5 


1 The MABL profiles correspondi ng to 
York sounding. (From Stage 1979) 
a. Temperature. 


b 


Mixing ratio. 
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data. 

The analytic solutions of horizontal heat and vapor 
transfer coefficients are shown in Fig. 2.5.2 
(Stage, 1983a, Fig. 10). The solutions are for the dry 
boundary layer. In this figure, the analytic solutions 
show that the horizontal heat (gi) and vapor (ga) transfer 
coefficients are 2.0 near the shore and eventually 
approach, zero when fetch (x) approaches inifinity. 

In order to evaluate the quality of the numerical 
solutions, they are compared with the analytical results. 
The numerical model has been run in both the wet 
(cloud formed) and the dry (no cloud formed) boundary 
cases, as shown in Fig. 2.5.3. The wet boundary layer 
case better resembles the real atmospheric situation* 
However, the result of the dry boundary layer case of the 
numerical solution is used to compare with the horizontal 
transfer coefficient solution which was analytically 
solved under the dry boundary layer condition. 

The so called numerical dry (cloud-off) case is an 
artificial situation, which was obtained by applying the 
dry (cloud-free) boundary equation. This equation 
forcibly suppresses the radiation and latent heat release 
terms when saturation is encountered. 

From the comparisons between the following two 
figures. Fig. 2.5.2 and Fig. 2. 5. 3, it is clearly seen that 
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Fig. 2. 5.2 The analytical solutions -for the dry 
MABL horizontal heat and vapor transfer 


coefficients. 
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the analytic solution, and the cloud free numerical 
solution of the horizontal heat and vapor transfer 
coefficients are identical at all fetches* This is 
theorectical ly accurate# Meanwhile, several conclusion 
can be made from these two figures, gi and g a as function 
of Log X. 

First, it can be inferred that once clouds form the 
entrainment rate increases due to radiation and latent 
heat release in the cloud layer, then the boundary layer 
begins dry out due to increased entrainment of dry air 
from above (Stage 1979, Fig# 5.2.7)# Therefore, once 
clouds form the boudary layer becomes drier than the 
cloud-off case at that same fetch. This will tend to 
increase the surface vapor fluxes and cause a large value 
of the horizontal vapor transfer coefficient for the 
cloud-on case. Also, the differences between the vapor 
transfer coefficients of the two cases increase following 
the boundary layer evolution. On the contrary, the 
horizontal heat transfer coefficient is larger for the 
cloud off case. It can be inferred too that once clouds 
form the boundary layer is more strongly heated by 
entrainment and is therefore warmer at any given fetch 
than the cloud-off case (Stage 1979, Fig. 5.2.8). This 
will tend to decrease the surface heat fluxes and cause a 
smaller value of the 'horizontal heat transfer 
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coefficient. 

Another noticeable point is that the analytical 
solutions start from zero fetch (x=0), i.e. Log X equals 
negative infinity. The numerical model sets up the 
initial condition of Z sz = 30m, with the corresponding 
value of Log X 2 -4.5. This assumption will be proven to 
be as good as a very thin boundary initial condition, eg: 
Z» z = 0.1 m, by using the modified computation of model 
initial condition in the next chapter, sec. 3.2. The 
advantage of this assumption is that it saves a lot of 
computation time during the model's initial stage* 

After all, it can be concluded that in the cloud-free 
regions the numerical solutions of horizontal transfer 
coefficients are equivalent to the analytical solutions 
derived by Stage (1983a). 

However, the biggest advantage of the numerical 
method is that it is good not only for the dry boundary 
layer, but it also provides a means of computing the 
horizontal transfer coefficients for the wet boundary 
layer evolution case, as has been shown in Fig. 2.5.3. 
Thus, this result can be used as a relatively accurate 
forecast estimates of surface heat and vapor fluxes from 
shore into the cloud-topped regions during cold-air 
outbreaks of marine boundary layer evolution. 

By comparing this section and last one, it is seen 
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that the horizontal transfer coefficient method for the 
boundary layer heat and vapor fluxes can be independently 
solved either numerically or analytically, as shown 
schematically in Fig. 2. 5. 4. 

In the next chapter, further studies of how the 
transfer coefficients respond to the variations of the 
fundamental parameters will be discussed in details. A 
modified technique of computing the heat and vapor fluxes 
into the cloud-topped boundary layer during cold air 
outbreaks is also developed in the following sections. 


2.6 Specification of The Model 

In the last section, the extended numerical version 
of Stage and Businger's boundary layer model (1981a, b) has 
been shown to be identical to the analytical solutions to 
the cloud edge. In next chapter, the model will be run to 
simulate evolution of the surface heat and vapor fluxes 
from shore into the cloud-topped regions. Those runs are 
done under wide ranges of change of the boundary layer 
fundamental parameters. The ranges of values have been 
chosen corresponding to the greatest uncertainty range of 
measured data likely to be encountered in the real 
atmosphere. Thus, the results of these runs could be 
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applied to the boundary layer evolution during most cold 
air outbreaks. 

This modified numerical model is formed under the 
assumption of linear shoreline sounding profiles in order 
to simplify the problem. Thus, the fundamental dependent 
variables $ m and q T of the boundary layer at the shore 
may be specified by the sounding values at the shoreline 
surface $ m * , g T i and the lapse rates Tl*; . 
Alternatively, the shoreline soundings could be specified 
in terms of linear and RH»,(i.e. in terms of Sxx, RH*, 
and the lapse rates T®i / Prh ) • Other initial values of 
variables are Z a = Z D = 0 , because we assume that Z B = Z. 
+ Z D = 0 , and Z B = Z. at the shore. 

Altogether, the full set of external parameters 
during a cold-air outbreak episode are qn, Ta-zHaT or 
flu, RHi, Toi r Trh which describe the shoreline sounding; 
wind speed (U); sea surface temperature (To); radiative 
sky temperature (T B kv), or cloud top radiation cooling 
rate (R*), and divergence (Di and D a ) * 

Thus, there are a total of 9 external parameters 
which are important to the evolution of the horizontal 
heat and vapor transfer coefficients during cold air 


outbreaks 
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2.7 Modification of Model Initialization 

The numerical model of simulating marine boundary 
layer evolution during cold outbreaks was first developed 
by Stage (1979). This model uses the assumptions of 
linear shoreline soundings as given in (2.4.1) and (2.4.2) 
and sets the initial values of fundamental variables as: 

Zbi = 0, 

$mx = 0mx, at X = 0 (2.7.1) 

<3tx = <3t i • 

9 m and q T are functions of Z B , which become independent of 
the initial conditions rapidly (Stage, 1983a). 

Stage (1979, 1983b) and Stage and Businger ( 1981a, b), 
were able to avoid spending large amounts of computation 
on the ' uninteresting and physically unrealistic initial 
growth of the boundary layer. They began their 

computation with 

Z b x = 30 m , 

$mx = Srn . B - - 0.001 > at x = 0 (2.7.2) 

q t i = qT , . 
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This technique was appropriate for their computation 
purpose* However, it can be seen from (2*3*1) and (2*3.2) 
that the net sensible (H) and latent (H c ) heat input to 
the boundary layer have different values for (2.7.1) and 
(2.7.2). For previous studies the difference is slight. 
Nevertheless, evolutions of H and He are the main topics 
of this paper, therefore a more accurate model 
initialization to estimate H and He is required here. 

The modified initialization method is based on the 
analytical solution (Stage, 1983a), and begins the 
numerical model computations at the fetch where the dry 
MABL was grown from the linear soundings of $ m and q T at 
the shore to Z B equal to 30 m. The net sensible (F 8IHB ) 
and latent heat (Flat) fluxes input are computed by the 
following steps. 

Z B i — 30 m 

B m Z = &m «. B-» 

(2.7.3) 

*Jti = <3t,b* 

Bit = Bx.u- = 0m. a* — ( L /Ci> ) qr , b + 

Assuming 0 X - 0* ( V q x s 0), from Stage (1983) 
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B » = + jfaA (Bx,u+ ~ $u) 


( 2 . 7.4 ) 


Then, F SCN b and Flat can be approximated by 9x . 


F ISMB = ( Bx 5" ( Bxx + )) Z, 


( 2 . 7 . 5 ) 


IT _ b / ^0 t 

- “cT F **~* < e.-e, A ' 


( 2 . 7 . 6 ) 


Because Bx - 9 , by using 


= ( 1 + ( 1/6 - 1 )g T1 ) B-. 


( 2 . 7 . 7 ) 


and bulk transfer formula, surface sensible heat flux, 
W'Pv'lo, can be computed. Then, the estimated time and 
fetch are obtained as: 


t = F, 


1INB / W ' S-v ' a , 


( 2 . 7 . 8 ) 


x = U t . 


( 2 . 7 . 9 ) 


Also, the potential temperature and mixing ratio of the 
boundary layer are estimated as: 


Flat 

•^ z b 


9 t = ( qi , b * 


+ qxi ) + 


f 


( 2 . 7 . 10 ) 



49 


Sm = 01 + - 7 =- qx . ( 2 . 7 . 11 ) 
'-P 

A further study of the horizontal tranfer coefficients, 
from (2.4.3) to (2.4.7), and the bulk transfer formulas, 
the initial values of 91 and g a can be derived. 


9i = 


CT/(ft©-&,)X 


( HW/C* ) 


C T , 




C T U( 6 10 * 1 )t 2 2.0 (2.7.12) 


ga = -T—Z 




( ) 




£77^0-?)* 


C^U (qo - 9 n )t ( L /Cp ) 2 2.0 (2.7.13) 


it is assumed that C T 5 C, and C T = C T 1 during the initial 
stage . 

This modified approximation of Fsinb and Flat during 
the initial stage of MABL evolution provides an explicit 
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set. of equations for model calculations and introduces 
only a very small errors. 

Meanwhile, it has also been proven that the initial 
analytical values of the horizontal transfer coefficients 
are approximated to be 2.0. 

It will be seen later in figures of Sec. 3*1 that the 
initial values of g* and g a of the numerical solution are 
almost identical to the analytical solutions derived here. 
This is another proof of the credibility of this modified 
numerical marine atmospheric boundary layer model. 


2.8 Modification of the Model Time Step 

This model uses a variable time step, DT . The time 
step depends on how rapidly the mixed layer evolves and is 
defined as : 


DT * 

where the default values of the parameters set up in this 
model are: dZmax = 20M; ddmax = 0.05; dqmax = 0.05xl0~ 3 
and fractional change of Z» = 0.05xZ». These constants 
limit the maximum changes in Z», Q m , q T and the fractional 
change of Z* which are allowed in any one time step of the 


rm* 

d9ntttx 

d9»/di 

' d9r/dt 


-Zfl • jmcr.frv 
' We+ fWiJ 


( 2 . 8 . 1 ) 


model run 
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In general, the boundary layer grows rapidly during 
the initial stage, so DT is small. Then, as time goes on 
the layer growth rate will decrease, and DT becomes larger 
and larger. If we use DT given by (2.8.1), trouble occurs 
for very long fetches, because DT is so large that the 
changes in Q m and q» in a single time step become too 
large as a fraction of the air-sea differences in $ m and 
g T . Meanwhile, the characteristics of the boundary layer 
evolution will be lost too.. 

In the extended version of this numerical model, 
additional limits on the size of the time step, DT, are 
imposed. In addition to satisfying (2.8.1), DT is limited 
by 



( 2 . 8.2 ) 


(2.8.3) 


These conditions guarantee that ff m0 ~$m and q T o - 3 t 
do not change by more than 10% in any one time step. 


2.9 Specification of Model Termination 

In addition, we put condition checks into the new 
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version of numerical model in order terminate the runs 
before physically unrealistic cases develop. These 
termination factors of the boundary layer evolution are 
shown and briefly explained in Table 2.9.1. 
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Table 2.9.1 

The terminating factors of the MABL model. 

I 

I 

Terminating factors | Physical reasons 

I 

I 

1. Z B > 5000 m When Z B becomes more than 5km, it 

is far beyond the normal BL height 
of 1-2 km. The well mixed boundary 
layer assumption can not be 
expected to hold for such large 
boundary layer depth. 

2. Zb < 30 m Z B has been unrealistically 

suppressed by divergence. 

3. Zc < 0 m This study does not include cases 

of stratus fogs. 

4. Tbkv < -273°C The radiative sky temperature 

must be above OK. 

5. Rb < 0 W/m a This study does not include cases 

of cloud top radiation warming 
during cold air outbreaks. 

6. X» < 0 Free entrainment cases are not 


included in this study 
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7. DTHV < 0 Cases of encroachment are not 

included in this study. 




CHAPTER 3 MODEL RESPONSES TO VARIATION OF EXTERNAL 
PARAMETERS 


3.1 The Responses of Horizontal Transfer Coefficients 

to The Variation of External Parameters 

Stage (1983b) has studied the evolution of basic 
thermodynamic variables of the marine boundary layer 
during cold air outbreaks. He pointed out how these 

factors were influenced by the variation of external 
parameters. This study will concentrate on how much these 
parameters influence the evolution of horizontal heat and 
vapor transfer coefficients. 

In this section, all the test runs are based on the 
New York sounding profile which has been studied by Stage 
and Businger (1981b), Chou and Atlas (1982), Stage 
( 1983a, b). This sounding is considered to be a typical air 
modification profile of cold-air outbreaks and is used to 
demonstrate the importance of various physical processes 
in determining boundary layer evolution. This sounding is 
shown in Table 3.1.1. 

In this section tests are run by varing each of the 
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Table 3.1.1 

Parameters of the basic test case ( based on 
sounding, Feburary 1979 ) 


Hater Surface 
Temperature 
Mixing Ratio 
Wind speed 
Divergence 

Radiative Sky Temperature 

Potential Temperature 
Sounding 

Total Hater Mixing Ratio 


To = 1 2 °C 

<3 O = Q B AT (To) = 

U =15 m/sec 

iMJ = 0 /sec 

T sky = -40 °C 
9 = 3.0 °C + 3.8 

<3t 


New York 


8 . 6g/kg 


C/km Z 


Sounding 


3 . 4g/kg-0 . 75g/kg/km Z 
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external parameters one at a time, and noting the 
resulting responses of the horizontal transfer 
coefficients • 


A. Wind Speed 

The test runs of wind speed varing from 1, 2, 4, 8, 
16 to 32 m/sec with the responses of g A and g a are shown 
in Fig *3. 1*1* 

As noted by Stage (1983a), wind speed has no effect 
on the distance required for cloud formation nor on the 
layer evolution before cloud forms* Low wind cases 
require a longer time to reach a given fetch, and cause 
smaller surface sensible and total water vapor fluxes, so 
that the radiation terms are dominant in the cloud regions 
for the low wind cases* In high wind cases, radiation is 
then relatively less dominant than surface fluxes. From 
(2*4*3), (2*4*4), (2*4*5) and the bulk transfer formulas: 


W ' io = CxU( Svo - 6 *. ), 


Wqx’lo = C„U( qxo - qx )• 


(3.1.1) 
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It can be shown that 

9i 50 ■ *° ^ v o - flv- ) dx , 

(3.1.2) 

/ * 

C <g ( <3tO - q T )dX, 
o 

and it is seen that g x and g a do not directly depend on 
wind speed in the dry case provided that C T and C q are not 
functions of wind speed (see discussions in Stage 1983a). 
However, from Fig. 3.1.1 it is seen that wind speed has 
little effect on g t and g a for the high wind speed 
(surface heat fluxes dominant) cases. For the low wind 
speed (radiation dominant) cases, although cloud top 
radiation cooling and surface heat fluxes input partly 
offset each other , g* and g a values are larger than those 
for higher wind speed. Also, the changes of g± and g a are 
much larger when wind speed are smaller, this suggetsts 
that wind speed has indirect effect on g x and g a . 

One important note is that with a large variation of 
wind speed ranging from 1 m/sec to 32 m/sec, the model 
predicted values of horizontal transfer coefficients gi 
and g a have less than 32% variation at any given fetch up 
to LOG X = 0 (i.e. fetch - 1900 km). In general g a is 
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less dependent on wind speed than gi, and more details are 
shown in Table 3.1.2* 


B. Divergence 

Two different types of divergence D t and D a are 
included. They were first defined in Stage (1983a). Di 
is the divergence confined in and resulting from the 
boundary layer. Changes of surface roughness and fluxes 
or the barocl inici ty in the boundary layer are all 
possible causes of Di, and are assumed to act only up to 
Z». Since this numerical model assumes a sufficiently 
convective boundary layer, the layer can be treated as a 
well mixed slab, Di can be assumed independent of height. 
D a is the synoptic scale divergence which is assumed to be 
uniform at all levels of interest, because the levels 
which are interested by our boundary layer studies are 
much less than the scale height for the synoptic weather 
patterns. 

The series of Di and D a runs in this section are done 
by setting cloud top radiation cooling rate (R») constant 
at 0.1 mK/sec (i.e. 128 W/m a ). 

First, let Di = -10x10“® to 10x10"® sec" 1 , in 
increments of 2x10"® sec" 1 , and D a = 0; i.e. from large 
convergence cases to divergence free and finally to large 
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Table 3.1.2 

The changes of g x and g 3 caused by a 32-fold change in 


wind 

speed 

. U : 1 

- 32 m/sec 

• 



LogX 

X 

Range of 

Range 

of 

Percent 

change in 



U 


9a 

9i 

9a 


(km) 

m/sec 

(9ii4gi) 

(ga+A9a ) 

(2igi/gi 

) ( 2A9a/9a ) 

-1.0 

193 

1 to 32 

1.49+0.11 

1.85+0.06 

14.0% 

6.7% 

in 

• 

0 

1 

611 

1 to 32 

1.1 9+0 . 19 

1.83+0. 15 

31.6% 

15.9% 

1 o 

1 • 

1 o 

1932 

8 to 32 

0.67+0.08 

1.63+0.08 

21.9% 

8.9% 
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divergence cases. The responses of gi and g a are shown in 
Fig. 3.1.2a. 

It is clearly shown from this figure that for any 
given fetch g x and g a differ by less than 55% of their 
range when Di is changed. The uncertainty of heat transfer 
coefficient g* is more than that of the mixing ratio 
transfer coefficient g a . It is also computed that the 
values of g* change less than 55% for the range of Ox 
values tested, and g a changes less than than 36%. More 
details are shown in Table 3.1.3a. 

The top of boundary layer is pushed upward by 
convergence and pulled downward by divergence, thus cloud 
can be formed at shorter fetches with larger convergence 
and suppressed with divergence. This is shown in 
Fig. 3.1.2b. For Di > 4xl0~® sec” 1 , no curves are shown 
in Fig. 3.1.2b, because a cloud can not be formed under 
such large divergence. 

Since the boundary layer becomes thinner with 
increased divergence, the mixed layer potential 
temperature 6m, and total water mixing ratio q T are 
higher, which has been shown by Stage (1983a, Fig. 8, and 
1983b Fig. 10). The effects of large divergence will 
increase the average $ m and g T of the boundary layer, 
although it is partially offset by the entrained cold and 
dry air from above cloud top , the surface heat and vapor 
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Tsbls 3*1* 3& 

The changes of g x and g a caused by Di type divergence. 
Di changes from -10x10*® to 10x10*® sec* 1 , in increments 


of 2x10“® 

sec* 1 , 

and 

o 

11 

II 

o 

• 




LogX 

X 

Range 

of 

Range 

of Percent 

change in 



Di 


9 1 

9 a 

9x 

ga 


(km) 

( 10“®sec* 1 ' ) 

(giiA9i) 

(9a+A9a) (2d9i/9i 

) (2Aga/?a) 

1 

• 

o 

193 

-10 to 

+ 10 

1.35+0.04 

1.81+0.06 

6.2% 

6.9% 

-0.5 

611 

- 2 to 

+ 10 

0.93+0. 12 

1.65+0.11 

24.7% 

13.8% 

o 

• 

o 

1932 

+ 4 to 

+ 6 

0.57+0.16 

1.35+0.24 

54.5% 

35.4% 




10 5 10 4 10 3 10 2 10 ' 

X 

Fig. 3.1«2b T 8 ky as function of 
divergence changes 
se 
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fluxes tend to decrease with the increased average 
potential temperature and mixing ratio of boundary layer. 
It is clearly seen on Fig. 3.1.2a that the horizontal 
transfer coefficients are smaller for larger divergences. 

Next, the effects of D a variation on g*, g a and T SK ir 
with constant R» at 0.1 mK/sec are shown in Fig. 3.1.3a 
and Fig. 3.1.3b respectively. During these test runs, D a 
= - 10x10”® to 10x10”® sec” 4 , in increments of 2x10”® 
sec" 4 , and D 4 = 0. 

The variations in gi and g a respond to changes in D a 
can be understood by using explanations used for D*. The 
ranges of variation in g* and g a due to the change of D a 
are less than 32% of the average values at any given 
fetch. The changes in g t are smaller than that of g a for 
the range of D a values tested. It is also noticed that 
the changes in g 4 and g a will tend to increase with 
increased fetch. Further details are discussed in Table 
3.1.3b as shown in the following page. 

In Fig. 3.1.3b, it is seen that the cloud is formed 
at a shorter fetch with stronger convergent motions. For 
D a larger than 4x10”® sec” 4 , no curves are shown in 
Fig. 3*1. 3b, because cloud can not be formed under 
suppression by such strong downward motion. 

From the conclusion above, we learn that basically 
the changes of Di and D a have similar effects on g t and 




X 

Fig. 3. 1.3a gi and g a as functions of X for D a type 


divergence changes. (D a : -10 to 10x10-° 

sec* 1 ) 
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Table 3.1.3b 

The changes of g*. and g a caused by D a type divergence. 
D a changes from -10x10“® to 10x10”® sec" 1 , in increments 
of 2x10"® sec" 1 and D* = 0. 


LogX x Range of 

Range 

Of 

Percent 

change in 

D a 

9i 

9a 

9i 

ga 

(km) (10"®sec" 1 ) 

(gi+Mi) 

( <3a+Aga ) 

(2dg»./g 

i) (2Ag a /5a) 


0 

• 

1 

193 

-10 

to 

+ 10 

1.31+0.03 

1.80+0.11 

4.3% 

12.6% 

1 

o 

• 

in 

611 

+ 2 

to 

+ 10 

0.75+0.10 

1 . 31 + 0 .21 

25.6% 

31.9% 

-0.37 

824 

+ 4 

to 

+ 8 

0.63+0.08 

1.30+0.16 

24.2% 

25.0% 
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g a . However, gi and g a do have a wider range of 
responses with respect to changes in the synoptic scale 
type divergence (D a ). Also, by comparing the values of g^ 
and g a at any given fetch for the same value of 0 t or Da? 
it is noticed that D a does have a stronger suppression on 
the vapor and heat fluxes input to the boundary layer as 
shown in Table 3.1.3a and 3.1.3b. 


C. Sea Surface Temperature (To) 

Test runs were done with To = 6, 8, 10, 12, 14, 16 
and 18 ° C, as shown in Fig. 3.1.4. The resulting 
variations of g t and g a due to water surface temperature 
(To) change are in a range of about 5 to 35% of its value 
at any given fetch. ga is much less sensitive to the 
change of To than g t is. More details are shown in 
Table 3.1.4. 

In Fig. 3;l.4, it is also seen that increasing To 
increases the surface heat fluxes, this tends to speed the 
growth , warming and moistening the boundary layer. 
Therefore, the boundary layer average temperature and 
mixing ratio will increase, which partially offset the 
increase of surface fluxes. This offset is larger for the 
heat fluxes, which can be seen in Fig. 3.1.4 that the 



Fig. 3. 1.4 g» and g a as -functions of X for 
various water surface temperature 
<T 0 * 6 to 18°C) . 





72 


Table 3.1.4 

The changes of g x and g a caused by a change in water 
surface temperature. To* T 0 : 6 to 18 °C, in increments of 
2°C. 

Log X x Range of Range of Percent change in 

To 9i 9a 9x 9a 

(km) (°C) (gi±A9i) (fa+dga) (2dgi/gi) (2Ag a /?a) 

-1.0 81-292 6 to 18 1.38±0.05 1.80+0.05 7.0% 5.3% 

-0.5 257-923 6 to 16 0.87+0.11 1.66+0.08 25.9% 9.6% 

0.0 1210-1932 8 to 12 0.60+0.11 1.55+0.02 36.8% 


2 . 1 % 




7.3 


* 

values of g* becomes smaller at large fetch for a warmer 
* 

water surface temperature. The cases with cooler water 
surface temperature have weaker surface fluxes and are 
more strongly affected by the radiation after cloud 
formed. The cloud top radiation cooling decrease the 
average boundary layer temperature and mixing ratio, hence 
the surface fluxes input become larger for the colder 
water surface temperature after cloud formation point. 
This phenomenon becomes more significant especially^ at 
very large fetch. It is also seen from Fig. 3.1.4 that the 
nondimenslonal fetch (X) starts off at a smaller value 
with larger water surface temperature. 

Another noticable feature is that sensible heat 
transfer coefficient gi decreases with increased To and 
partially offsets the sensible heat flux increase due to 
larger To* However, vapor transfer coefficient g a 
increases with increased To/ and thus enhances the vapor 
heat flux effect. 

D. Sky Temperature (T BK y) 

Changing radiative sky temperature from 0, -20, -40, 
-60 to -80 °C is one way to study the boundary layer 
radiation cooling effects on the cloud top. The test runs 
are shown in Fig. 3.1.5. First, we noticed that for this 
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2.0 


1.5 
1.0 
0.5 
0 

10~ 5 10' 4 10" 3 10~ 2 10 _1 10° 10 1 

X 

Fig. 3. 1.5 g t and g 3 aa functions of X for 
various radiative sky temperature' 
(TokvI 0 to -80 C) • 



i 
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range of change of T BK y the change of g* is between 5 to 
25% of its average value, and the range of change of g a is 
even smaller between 2 to 5% of its average value at any 
given fetch. Further details are presented in 
Table 3.1.5. 

Decreasing Tbky increases the cloud top radiation 
cooling rate, R» , so there is a larger entrainment rate 
and thicker cloud. The surface sensible and latent heat 
fluxes increase because of the decrease of average 
boundary layer potential temperature and total water 
mixing ratio, and they partially offset the radiation 
cooling effect. But, we can see from Fig. 3.1.5 that both 
heat and vapor horizontal transfer coefficients become 
larger with a smaller T B ky temperature. This means that 
the surface heat and vapor fluxes increase dominates the 
effect of increased radiation cooling and the effect of 
entraining cold and dry air from above cloud top into 
boundary layer. Therefore, it can be concluded that 
decreasing T BK v increases the cloud top radiation cooling, 
and it also speeds up the boundary layer growth. 

These runs are terminated when T B * becomes smaller 
than T BK y/ because this produces cloud top radiation 
warming, i.e. R B < 0, which is considered physically 
unrealistic in cold-air outbreaks. 
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Table 3.1.5 

The changes of g* and g a caused by Texv* Tsky! 0 to -80°C, 

0 

in increments of -20 C. 

Log X x Range of Range of Percent Change in 

Tbky 9i ga gi ga 

(km) ( a C) (9i+4gi) (ga+A9a ) ( 2^gi/gi ) ( 26g a /g 3 ) 

-1.0 193 -20 to -80 1.43+0.03 1.80+0.02 4.4% 1.8% 

-0.5 611 -20 to -80 1.03+0.07 1.72+0.04 13.8% 4.6% 


0.0 1932 -20 to -80 0.51+0.06 1.53+0.04 25.0% 


5.2% 
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E» Cloud Top Radiation Cooling (Rb) • 

Another way to specify the cloud top radiation 
cooling effect is to use constant value of R a through each 
run. R b of 0., 0.05/ 0.10/ 0.15/ 0.20/ 0.25/ 0.30 mK/sec 
were used in the test runs shown in Fig. 3.1.6. This 
corresponds to radiative heat fluxes from about 0 to 384 
W/m a in increments of 64 W/ro a . The ranges of variation of 
g» and g a are still quite small for the range of change of 
Rb being tested. Ranges of both horizontal heat and vapor 
transfer coefficients increase with increased fetch. Among 
the two coefficients/ g a is less sensitive to the 
variation of R a than is gi. Further details are shown in 
Table 3.1.6. 

Similar to the discussions of test runs for T BKyr it 
can be inferred that stronger cloud top radiation cooling 
can cause larger entrainment rate and speed up boundary 
layer growth. Therefore/ when R a = 0.25 and 0.30 mK/sec/ 
Zb exceeds 5km at very short fetch/ and the run are 
terminated. This is seen on Fig. 3.1.6. 

As has been shown in Fig. 3.5.2b and Fig. 3.5.3b/ it 
can be easily seen that the variations of radiative sky 
temperature (T B kv) were very large during boundary layer 
evolution while holding radiation cooling constant. So, 
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2.0 
1.5 
1.0 
0.5 
0 

10 -5 10~ 4 10~ 3 10 -2 10 _1 10° 10 1 

X 

Fig. 3. 1.6 gi and g 3 as functions of X for 
various cloud-top radiation cooling 
rate <R»i 0 to 256 W/m 3 ) . 
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T flbl 6 3 • 1 • 6 • 

The changes of gi and g a caused by R». R a : 0.0 to 0.30 
isK/sec, in increments of 0.05 mK/sec. 


LogX x Range of Range of Percent Change in 

RB 9i 9a 9i 9a 

(km) ( W/m a ) (gii49i) (9a+A9a ) ( 2Agi/gi ) ( 249a/9a ) 


-1.0 

193 

0 

to 

256 

1.39+0.06 

1.80+0.02 

8.0% 

1.8% 

-0.5 

611 

0 

to 

256 

0.95+0. 10 

1.72+0.05 

20.0% 

5.6% 

O 

• 

o 

1932 

0 

to 

128 

0.33+0.08 

1.46+0.06 

47.6% 

7.6% 
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holding T* KV constant is not a good assumption. 

However, physically holding cloud top radiation 
cooling (R») equal to constant is probably a better 
assumption than holding the radiative sky temperature 
constant as was done in the previous section. This 
concept can be easily explained. As the boundary layer 
evolves, the cloud top becomes higher and therefore 
colder. Meanwhile there is less water vapor above Z a and 
the air above Z a is colder, thus the radiative sky 
temperature also decreases. Since both the temperature at 
cloud top and the average radiative sky temperature above 
cloud top decrease, cloud top radiation cooling will 
remain relatively constant. 


F. Initial linear shoreline sounding ( 6 m ». , Vemi 
or ( RHx,T**h ) 

The initial linear shoreline sounding is a set of 
important external parameters for this model, we can 
either specify ( 9 ml , T9m} qn, Vmr ) or ( 0n, T* 0 i ; 
RH t , Tmm ) as described in sec. 2 . 6 . 

Fig. 3.1.7 shows model runs using linear soundings of 
potential temperature (0 m ) and total water mixing ratio 
(g T ) at shore. A series of 4 different sets of initial 


shoreline sounding runs has been done in order to see how 
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the evolution of boundary layer is affected by the 
shoreline sounding. They are: 

1. standard sounding except with $ ml increased by 2°C, 

2. standard sounding except with Ta* = 3.0*C/km, 

3. standard sounding except with g T x increased by lg/kg, 

4. standard sounding except with = 0 g/kg/km, 

where the standard sounding is the New York sounding as 
described in Table 3.1.1. 

In Fig. 3.1.7, it is seen that the initial points of 
all these runs start off at different locations; this can 
be inferred from (2.4.13), (2,4.14b), and (2.2.19) that 

nondimensional fetch is affected by the changes of initial 
sounding data of $ m i, qn, T$m , T«t. However, if the 
initial linear assumption of 6\ and RH is used, then the 
changes of RH*, f«i, and Phh can still affect and 

f w , and these two parameters are the key factors to the 
changes of x - and X. 

The changes of g t and g a due to initial sounding 
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parameters changes are in the range of 2% to 14% at any 
given fetch. Thus, it can be inferred that the horizontal 
heat and vapor transfer coefficients are not affected too 
much by the changes of initial soundings. More details 
are shown in Table 3.1.7. 

We may also use another set of linear shoreline 
sounding assumptions in which and RH are assumed to be 
linear at the shore, with surface liquid water potential 
temperature di i = 3.0°C at shore and the lapse rate = 

3.8 °C/km. This is useful because it allows us to check 
the effects of relative humidity by running the model with 
a series of different shoreline surface relative 
humidities, RHi = 0.2, 0.4, 0.6, and 0.8, with the lapse 
rate Frh = 0.0 g/kg/km. 

One thing noticed from Fig. 3.1.8 is that the 
variations of g* caused by the changes of RH* are much 
smaller than those of g a . From (2.4.8) to (2.4*16), g* is 
a function of X and is directly affected by RHj.. Since X 
becomes larger with a larger RHi from (2.4.12), so does g a 
and the changes of g a caused by RHi are significant. 
However, gi is not directly affected by RHi, the changes 
of RHi. cause variations of the rate of change of 5 v i, 
which will slightly change and g*. 
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Table 3.1.7 

The changes of g t and g a caused by the changes of 
shoreline surface initial sounding parameters. 

LogX x Range of Range of Percent change in 

initial g x g a , 9t g a 

sounding 

(km) parameters ( <?i±4gi ) ( 9a + 4g a ) ( 2Agi /gi ) ( 2ag a /g a ) 


-1.0 149-369 0.1: 1.39+0.02 1.74±0.02 2.9% 2.4% 

11.5-13.5C 
gn ! 

2.4-3.4g/kg 

-0.5 472-510 1.07+0.01 1.59+0.05 1.9% 6.5% 

( 10-*°C/km) 

1.925-3.0 

0.0 1494-1612 r*iT ’ 0.60 + 0.04 1.40 + 0.06 13.8% 8.1% 

( 10-»g/kg/km) 

-0.75-0.0 
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It is clearly seen that there are larger heat and 
vapor transfer coefficients at any given fetch when 
shoreline surface relative humidity is higher, this will 
tend to increase the surface heat and vapor fluxes. 
Meanwhile, the shoreline mixing ratio becomes larger with 
higher RHi, this tends to decrease the surface surface 
fluxes. These two effects partly offset each other, and 
the surface fluxes are less senstive to the changes of 
RH x . The variations of g a for the changes of RHi are much 
larger than the variations of g la Also, it is shown in 
Table 3.1.8 that the changes of both g x and g a increase 
with increased fetch. More details of the variations of 
gi and g a by the effects of shoreline surface relative 
humidity, RH 4 , are shown in Table 3.1.8* 


3.2 Summary of the Sensitivity Analysis 

The analysis done in the last section shows that gt 
and g a are not very sensitive to external parameters when 
those parameters change within the range likely to be 
encountered during cold-air outbreaks. It is therefore 
possible to obtain accurate approximation of heat and 
vapor fluxes even when the data contain large errors* 

Considering the New York sounding and a fetch of 600 
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Table 3.1.8 

The changes of gi and g a caused by RHi . RH*.: 0.2, 0.4, 
0.6, and 0.8, with x =3.0*C, T>i = 3.8®C/km, T*h = 0.0 
g/kg/km. 


LogX x Range of Range of Percent change in 

RHi gi g a 9i 9a 

(k») (9 i±A9i) (9a±A9a) (2A9 i/9i.) (249a/9a) 


-1.0 186-193 0.2 to 0.8 1.40+0.01 1.77+0.09 1.5% 10.6% 

-0.5 588-611 0.2 to 0.8 1.03+0.03 1.63+0.17 6.1% 20.4% 

0.0 1859-1933 0.2 to 0.8 0.55+0.03 1.38+0.18 11.3% 25.6% 




_as 


km (i.e. Log X s -0.5). Table 3.2.1 shows the range of 
variation of g x and g a when each parameter Is allowed to 
change by up to 50% from the New York sounding value. 

In this chapter. the major findings of my research 
have been presented. From the results of my work, it can 
be inferred that the surface heat and moist fluxes for 
either cloud-free or cloud-topped regions can always be 
reasonably estimated by knowing the analytical values of 


gw ga at 

any fetch, the 

shoreline 

temperature 

lapse 

rate. and 

the surface 

temperature 

and 

mixing 

ratio 

differences 

between air 

over land 

and 

sea • 

This 

approximation should be 

reasonable 

accurate for 

all 


different soundings. 

3.3 The Effect of External Parameters on the Horizontal 
Transfer Coefficients 

In this section we look at the first (direct) and 
second (indirect) order effects of all external parameters 
on the horizontal transfer coefficients. 

The independent variable of this marine atmosphere 
boundary layer (MABL) model is either time (t). or 
distance of fetch following the mean motion of air in the 
boundary layer (x). The relationship between t and x is: 
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Table 3.2.1. 

The uncertainties of g* and g» caused by the errors of the 
external parameters at fetch equal to 600 km (The 
measurement errors are allowed to be up to 50%). 


external data 

tolerance 

range of 

resulting 

errors 

(w/ standard 

of data 

parameters 

(fetch = 

600km ) 

sounding 



(gi+Agi) 

(ga+oga ) 

values 



( 24gi/gi) 

( 2 Ag a /5a ) 


Initial sounding: 


0 e i (*C) 11.5, 

q-ri (g/kg) 3.4, 
T* a (°C/km) 3.8, 
r, T (JgsL) -0.75, 


( changes of ) 


sounding 

data 


11.5 

to 

13.5 


3.4 

to 

4.4 1.07+0.01 

1.59+0.05 

3.0 

to 

( 1.9%) 

3.8 

( 6.3%) 


-0.75 to 0.0 
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Table 3 . 2 . 1 . (continued ) 


external data 
(w/ standard 
sounding 
values 

tolerance 
of data 

range of 

parameters 

resulting 
(fetch = 
(gi±A9i ) 

( 2agi/gi ) 

errors 
600km ) 
(9a±A9a ) 

( 2ag a /ga ) 

u 

(ro/sec) 16, 

50% 

8 

to 

24 

1 .01±0.02 

1.73+0.02 







( 3.9%) 

( 2.3%) 

Di 

( 10“®/sec ) 4 , 

50% 

2 

to 

6 

0.96+0.06 

1.58+0.06 







( 12.5%) 

( 7.6%) 

D a 

( 10~®/sec ) 4 , 

50% 

2 

to 

6 

0.84+0.03 

1.41+0.10 







( 7.1%) 

(13.5%) 

To 

(°C) 12, 


8 

to 

16 

1.03+0.11 

1.72+0.03 







(21.4%) 

( 3.5%) 

Tsk v ( ° C ) -40, 


-20 

to 

-60 

1.00+0.04 

1.72+0.02 







( 8.0%) 

( 2.9%) 

Re 

(mK/sec) 0.1, 

50% 

.05 

to 

. 15 

1.01+0.05 

1.72+0.03 







( 9.9%) 

( 3.5%) 

RH i 

(%) 40, 

50% 

20 

to 

60 

1.03+0.03 

1.55+0.08 







( 6.1%) 

( 10.7%) 




9.1 


x = U t , 


(3.3.1) 


where 


1) is the mean wind speed in MABL. 

All the other dependent variables of the MABL are 
functions of either time or fetch. We will now look at the 
model evolution of these basic parameters of the boundary 
layer^ e.g. $ m , q» • 

First, the potential temperature change per unit time 
can be derived from (2.3.13)/ (2.3.14)/ and (2.3.16). 

-gl 2 - =-^( W'0. '|o + R c - R> + W.A*. ) (3.3.2) 

where 

w*0. 4o = CxU ( 9 mo - 9m ) (3.3.3) 


In another expression/ the temperature change per unit 
travel can be expressed as: 


_ J_ dPa [_ , w'&a' o Weggg K* _ Ra \ 

dx u dt ’ Zb U u u “ U ; 
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Then, from (3*3.3), (3.3.4) can be written as: 



~~ ( Ct ( $ mO~ 8m ) ♦ 
*a 


u 



) . 


(3.3.5) 


From (3.3.2) to (3.3.5), it is proven that wind 
speed (U) has first order effect on the changes of 
temperature per unit time and travel; surface heat fluxes 
per unit time; radiation warming, cooling and entrained 
heat fluxes per unit travel. However, the wind speed has 
only second order effect on surface heat fluxes per unit 
trave 1 . 

Similarly, the rate of change of total water mixing 
ratio per unit time can be derived from (2.3.12), and 
(2.3.15) . 

-gj 1 - = ( vFi^’lo + W.Agx ), (3.3.6) 


where 


Wqx'lo = C,U( qo - qx ) , 


(3.3.7) 


and then the mixing ratio change per unit travel can be 


expressed as 
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d<\r 

dx 


UW S % <C KV* > * 


W^i- ) 
u ' 


(3.3.8) 


A similar conclusion could be obtained using (3.3.7) 
through (3.3.10). The wind speed has first order effect 
on the change of total water mixing ratio per unit time 
and travel, surface water vapor fluxes per unit time, and 
entrained vapor fluxes per unit travel. However, it has 
second order effect on the surface water vapor fluxes per 
unit travel. 

From, (2.4.3), (2.4.4), and (2.4.5), g* and g a can 
also be expressed as: 




fCp 

Ci 



6*o- (9a ) 


dx , 


(3.3.9) 


9a 


SL 

C*X<V<fT,) 



(q»-q-r 


dx , 


(3.3.10) 


As is concluded that wind speed does not directly 
affect the surface heat and vapor fluxes per unit travel. 
It is also seen from (3.3.9) and (3.3.10) that gi and g a 
are not directly affected by wind speed. However, wind 


speed does indirectly affect g t and g a , 


because the rates 
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of change of 9 m and g T are affected by wind speed (from 
(3.3.2) and (3.3.6)). This has also been identified in 
Fig. 3.1.1, which shows that when wind speed is large the 
changes of g 4 and ga are small with respect to the 
variations of wind speed; and when wind speed is small, 
the changes of g & and g a are much larger. There is no 
linear relationship between wind speed and horizontal 
transfer coefficients. This also tells that wind speed 
has second order effect on the horizontal transfer 
coef f icients . 

In a more generalized conclusion, from (3.3.2) and 
(3.3.6), it can be inferred that all the external 
parameters have a second order effect on the surface heat 
and vapor fluxes per unit travel, i.e. the horizontal 
transfer coefficients. 

In the next section, the strengths and weaknesses of 
this numerical boundary layer method will be discussed, 
and some further research direction in this area will be 
discussed in the next chpater. 


3.4 Model Strengths and Weaknesses 

In Chapter 2 and 3, we have compared the results of 
analytical and numerical solutions of the marine 


atmosphere boundary layer evolution and studied 


the 
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evolution of the horizontal transfer coefficients with 
respect to the boundary layer external parameter 
variation. In conclusion, the strengths and weaknesses of 
the modified numerical model will be discussed in this 
section . 

The analytical solution derived by Stage (1983a) is 
valid only for the dry boundary layer, cloud-free cases. 
While the numerical integration method presented in this 
paper not only provides an identical solution to the 
analytic solution of dry boundary layer case (sec. 2.5), 
it also presents a method to estimate the solution for 
the wet, cloud formed, boundary layer condition. 

Hence, the most significant point of the 
theoretically accurate numerical method is that it 
provides an easy way of computing the MABL heat and vapor 
fluxes during cold-air outbreaks. It is good for both 
before and after cloud formation during boundary layer 
evolution, and it is much more applicable to the real 
atmosphere stratus topped boundary layer conditions than 
is the analytic model. 

Another noticeable strength of this numerical model 
is that it can predict the estimated boundary layer 
evolution from conveniently available data sources, such 
as the shoreline sounding parameters ((Bai,!)*, qn/T«T) , 
wind speed (U), divergences (Di, Da)), and satellite 
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measurable parameters (water surface temperature (To)» and 
radiative sky temperature (T« K v))* 

The changes of the predicted values of horizontal 
heat and vapor transfer coefficients are quite small with 
respect to the wide range changes of measured external 
boundary layer parameters, as has been shown in Section 
3.5. 

Therefore, by using this numerical method, we can 
afford to have large measurement errors of the sounding or 
satellite data, and still come out with reasonably 
accurate estimates of the boundary layer heat and vapor 
fluxes input from the warm water surface. 

On the other hand, because this model is based on the 
idealized initial condition for analytic solution, there 
are also some constraints on this numerical method. As 
mentioned in Chapter 2, the idealized initial assumptions 
of horizontal homogeneous state in the shoreline 
direction, linear boundary layer profiles often do not 
exist in real atmospheric boundary layers. 

Also, if the contrast in temperature between the 
water surface and the air is not large enough, buoyant 
production is no longer dominant, and a more complicated 
MABL model which includes the shear production will be 
needed. The horizontal homogeneous assumption in the 
shoreline direction is also quite severe. It 
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oversimplifies the boundary layer evolution in the 
shoreline direction, and the computation of surface heat 
and vapor fluxes per unit area from the warm water 
surface . 

Furthermore, this model is not complete enough to 
include the conditions of free entrainment on the top of 
the boundary layer, radiation warming at cloud top, and 
radiation cooling at the cloud base. As mentioned in 
Section 3.4, this numerical model will be terminated once 
these additional situations above are encountered. 

For further studies in this area of MABL evolution 
during cold air outbreak episodes, all the critical 
conditions discussed above must be included in this model 
also. Therefore, further development of a more 
generalized complete numerical model to closely simulate 
the real atmospheric situations need to be done in future 


endeavors 



CHAPTER 4 


Conclusions and Outlook 


Stage and Businger ( 1981a, b) first proposed the MABL 
model which simulates the growth and evolution of a cloud- 
topped boundary layer during cold-air outbreaks. This 
thesis using the modified version of Stage and Businger's 
model numerically studies the surface heat and vapor 
fluxes into the boundary layer. Several significant 
results have been found through the study of this paper. 

Chou and Atlas (1982) first proposed and proved the 
idea of horizontal transfer coefficients for the cloud- 
free region. This is a new kind of parameterization for 
the convectively unstable boundary layer case, which 
depends on the horizontal temperature and mixing ratio 
differences across the coast, instead of the vertical 
differences as in the classical bulk transfer formulas. 
Through the present work this idea has been extended from 
cloud-free regions into cloud-topped regions. The 
horizontal heat (g t ) and vapor (g a ) transfer coefficients 
numerically estimated by the modified MABL model have been 
proven to be identical to the analytical solutions by 
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Stage (1983a) for the dry (cloud-free) case, and also can 
be applied to the wet (cloud-topped) case. Thus, the 
numerical MABL model simulates the marine boundary layer 
surface fluxes more realistically than the analytical 
method does when stratus clouds are present in the upper 
boundary layer. In addition, from my work it is found 
that the changes of g* and g a are not very sensitive to 
the changes of those sounding memasurable external 
parameters . 

Nevertheless, the most significant conclusion can be 
inferred that the estimates of gi and g a at any given 
fetch should always be reasonably accurate for any given 
soundings by using the analytical (dry boundary layer) 
solutions figure of gi and g a (Fig. 2.5.2, or Fig. 2.5.3). 
The values of fetch can be computed from the shoreline 
sounding temperature lapse rate, and the surface 
temperature differences between air at shore and sea. 
Hence, the surface heat and vapor fluxes at any given 
fetch and the total fluxes from shore to any value of 
fetch in the boundary layer, cloud-free or cloud-topped, 
are easily obtained. 

Since only a few initial shoreline sounding data are 
required for the run of simulating the boundary layer 
evolution, it is also recommended to actually run this 
MABL model numerically if more accurate results are 
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desired. Therefore/ the actual values of g* and g a and 
total fluxes can be obtained from the numerical model 
results, which are very good estimates even with some 
large sounding measurement errors as discussed early in 
Chapter 3. 

The boundary layer is an important link of energy 
budgets between the oceanic and atmospheric circulations. 
Further work in this area should include a more detailed 
evaluation of distribution of the energy source from the 
water (i.e. the buoyancy production of Turbulent Kinetic 
Energy); how much of this energy is dissipated by the 
eddies of the boundary layer; how much is used to drive 
the entrainment process; and how much is carried up beyond 
the boundary layer to serve as an energy input for the 
atmosphere circulation. Atlas and Chou's (1983) studies 
of the feedback effects from ocean to the low level 
mesoscale circulation in the Northern Atlantic winter 
suggested that the atmospheric circulations are greatly 
affected by the boundary layer surface heat and vapor 
fluxes from water. Further studies in this area will be 
very useful for the energetics of cyclogenesis during 
winter time storms over warm water surfaces, such as 
Northern Atlantic, Northern Pacific ocean and East China 
sea winter low pressure system regions. 

A more generalized MABL model which should include 
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not only the buoyancy production, but also the shear 
production. This model could simulate the boundary layer 
evolution for any well mixed layer, i.e. the buoyancy 
dominant unstable layer, or the near neutral layer as long 
as the shear term is strong enough to keep the layer well 
mixed. 

Overall, the most important topic for future research 
should be emphasized on how the oceanic heat and vapor 
fluxes contribute to the energy and vorticity of 
atmosphere circulation by passing through boundary layer 
during the cold-air outbreaks episodes. A more complete 
numerical model including the Turbulent Kinetic Energy 
equation and the Turbulent Vorcity Budget equation might 
be able to offer some clues to solve this problems of 


future endeavors 
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